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^ ■ Abstract 

This is a review of the dynamics of wave propagation through a disordered 

iV-mode waveguide in the localized regime. The basic quantities considered 

are the Wigner-Smith and single-mode delay times, plus the time-dependent 

power spectrum of a reflected pulse. The long-time dynamics is dominated 

by resonant transmission over length scales much larger than the localization 

length. The corresponding distribution of the Wigner-Smith delay times is 

i-rt ! the Laguerre ensemble of random-matrix theory. In the power spectrum the 

resonances show up as a t~ 2 tail after N 2 scattering times. In the distribu- 

q ■ tion of single-mode delay times the resonances introduce a dynamic coherent 

backscattering effect, that provides a way to distinguish localization from ab- 

C\| . sorption. 
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O 
O ! I. INTRODUCTION 

o 

Light localization, one of the two central themes of this meeting, has its roots in electron 
localization. Much of the theory was developed first for electrical conduction in metals at 
low temperatures, and then adapted to propagation of electromagnetic radiation through 
disordered dielectric media |l],|[. Low-temperature conduction translates into propagation 
that is monochromatic in the frequency domain, hence static in the time domain. 

This historical reason may explain in part why much of the literature on localization 
of light deals exclusively with static properties. Of course one can think of other reasons, 
such as that a laser is a highly monochromatic light source. It is not accidental that one of 
the earliest papers on wave localization in the time domain || appeared in the context of 
seismology, where the natural wave source (an earthquake or explosion) is more appropriately 
described by a delta function in time than a delta function in frequency. 

Our own interest in the dynamics of localization came from its potential as a diagnostic 
tool. The signature of static localization, an exponential decay of the transmitted intensity 
with distance, is not unique, since absorption gives an exponential decay as well M. This is at 
the origin of the difficulties surrounding an unambiguous demonstration of three-dimensional 







FIG. 1. The top diagram shows the quasi-one-dimensional geometry considered in this review. 
The waveguide contains a region of length L (dotted) with randomly located scatterers that reflects 
a wave incident from one end (arrows). The number of propagating modes N may be arbitrarily 
large. The one-dimensional case N = 1 is equivalent to the layered geometry shown in the bottom 
diagram. Each of the parallel layers is homogeneous but differs from the others by a random 
variation in composition and/or thickness. 

localization of light ||. The dynamics of localization and absorption are, however, entirely 
different. One such dynamical signature of localization |J is reviewed in this lecture. 

Localization is a non-perturbative phenomenon and this severely complicates the theo- 
retical problem. In two- and three-dimensional geometries (thin films or bulk materials) not 
even the static case has been solved completely M. The situation is more favorable in a 
one-dimensional waveguide geometry, where a complete solution of static localization exists 
]7|||. The introduction of dynamical aspects into the problem is a further complication, and 
we will therefore restrict ourselves to the waveguide geometry (see Fig. 1). The number N 
of propagating modes in the waveguide may be arbitrarily large, so that the geometry is 
more appropriately called quasi-one- dimensional. (The strictly one-dimensional case N — 1 
is equivalent to a layered material.) 

The basic dynamical quantity that we will consider is the auto-correlator of the time- 
dependent wave amplitude u(i), 

/oo 
dt'e- wt u{t)u{t + t'). (1) 

-oo 

If the incident wave is a pulse in time, then the transmitted or reflected wave consists of 
rapid fluctuations with a slowly varying envelope (see Fig. 2). The correlator a u (t) selects the 
frequency component u of the rapid fluctuations. The remaining ^-dependence is governed 
by the propagation time through the waveguide. 

If the incident wave is not a pulse in time but a narrow band in frequency, then it is 
more convenient to study the frequency correlator 

a w {5u) = / dte i&u)t a w (t) = u*{uj)u(uj + 8u). (2) 
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FIG. 2. Computer simulation of an acoustic plane wave pulse reflected by a randomly layered 
medium. The medium is a model for the subsurface of the Earth, with a sound velocity that 
depends only on the depth. The figure shows the reflected wave amplitude as a function of time 
(arbitrary units). The incident pulse strikes the surface at time zero. From Ref. M. 
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FIG. 3. Frequency dependence of the phase (modulo 2tt) of microwave radiation transmitted 
through a disordered waveguide. The waveguide consists of a 1 m long, 7.6 cm diameter copper 
tube containing randomly positioned polystyrene spheres (1.27 cm diameter, 0.52% volume filling 
fraction). Wire antennas are used as the emitter and detector at the two ends of the tube. From 
Ref. 0. 

The Fourier transformed wave amplitude u(u>) = J dt e llJjt u(t) = I l ' 2 e^ is complex, contain- 
ing the real intensity I{ui) and phase 4>(ui). Most of the dynamical information is contained 
in the phase factor, which winds around the unit circle at a speed d<p/duj determined by the 
propagation time (see Fig. 3). 

The correlator a depends sensitively on the random locations of the scatterers in the 
waveguide, that give rise to the localization. This calls for a statistical treatment, in which 
we consider the probability distribution of a in an ensemble of waveguides with different 
scatterer configurations. The method of random-matrix theory has proven to be very effec- 
tive at obtaining statistical distributions for static scattering properties ||. The extension 
to dynamical properties reviewed here is equally effective for studies of the reflected wave. 
The time dependence of the transmitted wave is more problematic, for reasons that we will 
discuss. 



II. LOW-FREQUENCY DYNAMICS 

The low-frequency regime is most relevant for optical and microwave experiments 
p yToyTTH , where one usually works with an incident beam that has a narrow frequency 
bandwidth relative to the inverse propagation time through the system. We assume that 
the length L of the waveguide is long compared to the (static) localization length £ = Nl, 



which is equal to the product of the number of propagating modes N and the mean free 
path I. The reflected wave amplitudes r mn in mode m (for unit incident wave amplitude in 
mode n) are contained in an N x N reflection matrix r. This matrix is unitary, provided we 
can disregard absorption in the waveguide. It is also symmetric, because of reciprocity. (We 
do not consider the case that time-reversal symmetry is broken by some magneto-optical 
effect.) 

The correlator 

C u (6u)=r\u)r(u + 6u) (3) 

is the product of two unitary matrices, so it is also unitary. Its eigenvalues exp(z</>„), n = 
1, 2, ... N, contain the phase shifts <fi n . Since <p n = for all n if Sou = 0, the relevant 
dynamical quantity at low frequencies is the limit 

r n = lim — -, (4) 

which has the dimension of a time. It is known as the Wigner-Smith delay time, after 
the authors who first studied it in the context of nuclear scattering fT2yT3|j . The r n 's may 
equivalently be defined as the eigenvalues of the Hermitian time-delay matrix Q, 

Q(oo) = -ir+— = f/ t diag(r 1 , r 2 , . . . r N )U. (5) 

Experiments typically measure not the product of matrices, as in Eq. @, but the product 
of amplitudes, as in Eq. (^|). The amplitude measured within a single speckle (or coherence 
area) corresponds to a single matrix element. The typical observable is therefore not the 
Wigner-Smith delay time but a different dynamical quantity called the single- channel (or 



single-mode) delay time |10| , pT 



T ™ rV m 5u\r mn ^)\ 2 ' { ' 

If we decompose the complex reflection amplitude into intensity and phase, r mn = J 1//2 e l<?i , 
then the single-channel delay time is the phase derivative, r mn = dip/duo = <\j ' . Since the 
reflection matrix r(u + 5uo) has for small 5a; the expansion 



(OJ 



+ 5u) = Y^ U km U kn {\ + ir k 5u), (7) 



we can write the single-channel delay time as a linear combination of the Wigner-Smith 
times, 

Trnn = <\> = Re— -, A k = ^ T*U im U in . (8) 

A i 

We will consider separately the probability distribution of these two dynamical quantities, 
following Refs. [|,|TJ. 



A. Wigner-Smith delay time 

There is a close relationship between dynamic scattering problems without absorption 
and static problems with absorption \T<^ . Physically, this relationship is based on the notion 



that absorption acts as a "counter" for the delay time of a wave packet |16| . Mathematically, 
it is based on the analyticity of the scattering matrix in the upper half of the complex plane. 
Absorption with a spatially uniform rate l/r a is equivalent to a shift in frequency by an 
imaginary amount 5u> = i/2r a .[] If we denote the reflection matrix with absorption by 
r(w, r a ), then r(u, r a ) = r{ui + i/2r a ). For weak absorption we can expand 



r(u + i/2r a ) « r(w) + ——r(uj) = r{u) 



1 " o-0(«) 

2r a 



(9) 



As before, we have assumed that transmission can be neglected so that r is unitary and Q 
is Hermitian. Eq. (|9|) implies that the matrix product rr^ for weak absorption is related to 
the time-delay matrix Q by a unitary transformation []H |, 



r(w,r a )r t (a;,r a ) =r(w) 



-QH 



r'(o;) 



(10) 



The eigenvalues R\, R2, ■ ■ ■ Rn of rr^ in an absorbing medium are real numbers between 
and 1, called the reflection eigenvalues. Because a unitary transformation leaves the eigen- 
values unchanged, one has R n = 1 — r n JT & . This relationship between reflection eigenvalues 
and Wigner-Smith delay times is useful because the effects of absorption have received 
more attention in the literature than dynamic effects. In particular, the case of a single- 
mode disordered waveguide with absorption was solved as early as 1959, in the course of a 
radio-engineering problem [IT] . The multi-mode case was solved more recently [|18| , |T9| . The 



distribution is given by the Laguerre ensemble, after a transformation of variables from R n 
to A n = R„(l — R n )~ '■ 



P({^n}) cx n |Ai - Xjfl[exp[-(aT B /T*)(j3N + 2 - (3)X k \. 



(11) 



i<] 



Here r s is the scattering time of the disorder and a is a numerical coefficient of order unity.0 
The symmetry index (5 = 1 in the presence of time-reversal symmetry. (The case /3 = 2 of 



1 To see this, note that absorption is represented by a positive imaginary part of the dielectric 
constant e = 1 + j/wr a (for ujt & S> 1). Since e is multiplied by uj 2 in the wave equation, a small 
imaginary increment u> — > uj + i8u is equivalent to absorption with rate 26uj. In the presence of 
a fluctuating real part of e, an imaginary shift in frequency will lead to a spatially fluctuating 
absorption rate, but this is statistically equivalent to homogeneous absorption with an increased 
scattering rate. 

2 The coefficient a depends weakly on N and on the dimensionality of the scattering: a = 2 for 
N = 1; for iV — > 00 it increases to 7r 2 /4 or 8/3 depending on whether the scattering is two or 
three-dimensional. The mean free path I, that we will encounter later on, is defined as I = o.'ct s , 



broken time- reversal symmetry is rarely realized in optics.) The eigenvalue density is given 
by a sum over Laguerre polynomials, hence the name "Laguerre ensemble" |M| . 

The relationship between the reflection eigenvalues for weak absorption and the Wigner- 
Smith delay times implies that the r n 's are distributed according to Eq. (|TTD if one substitutes 
An/r a — > l/r n (since X n — > (1 — i?„) _1 for weak absorption). In terms of the rates \i n = l/r„. 



one has [14 



P(M) oc 1] \fk - )"/ II ex P[-7(^ + 2 - /3)H- (12) 

«<j k 

We have abbreviated 7 = «r s . For iV = 1 it is a simple /3- independent exponential distribu- 



tion |2T| p3|, or in terms of the original variable t, 

P(t) = 27r- 2 exp(-2 7 /r). (13) 

The slow t~ 2 decay gives a logarithmically diverging mean delay time. The finite localiza- 
tion length £ is not sufficient to constrain the delay time, because of resonant transmission. 
Resonant states may penetrate arbitrarily far into the waveguide, and although these states 
are rare, they dominate the mean (and higher moments) of the delay time. The divergence 
is cut off for any finite length L of the waveguide. Still, as long as L > ^, the resonant 
states cause large sample-to-sample fluctuations of the delay times. These large fluctuations 
drastically modify the distribution of the single-channel delay time, as we will discuss next. 

B. Single- channel delay time 

In view of the relation (|8|), we can compute the distribution of the single-channel delay 
time <f)' from that of the Wigner-Smith delay times, if we also know the distribution of the 
matrix of eigenvectors U. For a disordered medium it is a good approximation to assume that 
U is uniformly distributed in the unitary group, independent of the r n 's. The distribution 
P (</)') may be calculated analytically in the regime N ^> 1, which is experimentally relevant 



(N ~ 100 in the microwave experiments of Ref. ||11| ). 

In the large- iV limit the matrix elements U mn become independent complex Gaussian 
random numbers, with zero mean and variance (\U mn \ 2 ) = 1/N. Since Eq. (RI) contains the 
elements U im and U in , we should distinguish between n = m and n ^ m. Let us discuss 
first the case n^mof different incident and detected modes. The average over the £/ in 's 
amounts to doing a set of Gaussian integrations, with the result 

P{#) = (\{B 2 - B\){B 2 + <P' 2 - 2iWr 3/2 >- (14) 



with a' = 2 for N = 1 and a' — ► tt/2 or 4/3, respectively, for N — > 00 in two or three dimensions. 
(The wave velocity is denoted by c.) Finally, the diffusion coefficient D = c 2 T s /d with d = 1 for 
N = 1 and d — ► 2 or 3 for N — > 00. The dimensionality that determines these coefficients is a 
property of the scattering. It is distinct from the dimensionality of the geometry. For example, a 
waveguide geometry (length much greater than width) is one-dimensional, but it may have d = 3 
(as in the experiments of Ref. pH ) or d = 2 (as in the computer simulations of Ref. §). 
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FIG. 4. Distribution of the single-channel delay time (j)' in the diffusive regime (top panel) and 
localized regime (bottom panel). The results of numerical simulations (data points) are compared 
to the predictions ( |14| ) (solid curve) and ( |T5| ) (dashed). These are results for different incident and 
detected modes n ^ m. From Ref. Q. 

The average (• • ■) is over the two spectral moments B\ and B 2 , defined by Bk = J2i T t\Uim\ 2 - 
The joint distribution P(B 1 , B 2 ), needed to perform the average, has a rather complicated 
form, for which we refer to Ref. M. 

The result ( |i~4"D applies to the localized regime L ^> £. In the diffusive regime I <C L <C £ 
one has instead fllTLEll 



P{#) = (Q/20)[Q + W/p - l) 2 ]- 3/2 - 



(15) 



The constants are given by Q ~ L/l and 0' ~ L/c up to numerical coefficients of order 
unity. Comparison of Eqs. flT4D and (O) shows that the two distributions would be identical 
if statistical fluctuations in B\ and B 2 could be ignored. However, as a consequence of the 
large fluctuations of the Wigner-Smith delay times in the localized regime, the distribution 
P(Bi, B 2 ) is very broad and fluctuations have a substantial effect. 

This is illustrated in Fig. 4, where we compare P(<j)') in the two regimes. The data points 
are obtained from a numerical solution of the wave equation on a two-dimensional lattice, 
in a waveguide geometry with N = 50 propagating modes. They agree very well with the 
analytical curves. The distribution fllB]) in the diffusive regime decays oc |0'| -3 , so that the 
mean delay time is finite (equal to 0'). The distribution in the localized regime decays more 



slowly, 



oc 



The resulting logarithmic divergence of the mean delay time is cut off in 
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FIG. 5. Same as the previous figure, but now comparing the case n/mof different incident 

and detected modes (solid circles) with the equal- mode case n = m (open circles). A coherent 

backscattering effect appears, but only in the localized regime. From Ref. |6|. 



the simulations by the finiteness of the waveguide length. 

Notice that, although the most probable value of the single-channel delay time is pos- 
itive, the tail of the distribution extends both to positive and negative values of <fi' . This 
is in contrast to the Wigner-Smith delay time r n , which takes on only positive values. The 
adjective "delay" in the name single-channel delay time should therefore not be taken lit- 
erally. The difficulties in identifying the phase derivative with the duration of a scattering 
process have been emphasized by Biittiker |25| . 

We now turn to the case n = m oi equal-mode excitation and detection. An interesting 
effect of coherent backscattering appears in the localized regime, as shown in Fig. 5. The 
maximal value of P{4>') for n = m is larger than for n ^ m by a factor close to \[2. (The 
precise value in the limit A" — > oo is v2 x ^^ ■) In the diffusive regime, however, there is 
no difference in the distributions of the single- channel delay time for n = m and n ^ m. 

Coherent backscattering in the original sense is a static scattering property [p6| , p7[ . The 
distribution P(I) of the reflected intensity differs if the detected mode is the same as the 
incident mode or not. The difference amounts to a rescaling of the distribution by a factor 
of two, 



P(I) 



-NI 



Ne 

l Ne -NI/2 [in 



if n 7^ m 



m 



(16) 



so that the mean reflected intensity I becomes twice as large near the angle of incidence. 
It doesn't matter for this static coherent backscattering effect whether L is large or small 
compared to £. The dynamic coherent backscattering effect, in contrast, requires localization 
for its existence, appearing only if L > £. This is the dynamical signature of localization 
mentioned in the introduction. 



C. Transmission 

Experiments on the delay-time distribution have so far only been carried out in trans- 
mission, not yet in reflection. The distribution fll5|) in the diffusive regime applies both to 
transmission and to reflection, only the constants Q and 0' differ |24|1 . (In transmission, Q is 



of order unity while <p' ~ L /lc.) Good agreement between theory and experiment has been 



obtained both with microwaves |TT| and with light |4J|. The microwave data is reproduced 
in Fig. 6. Absorption can not be neglected in this experiment (L exceeds the absorption 
length / a by a factor 2.5), but this can be accounted for simply by a change in Q and </>'. The 
localization length is larger than L by a factor of 5, so that the system is well in the diffusive 
regime. It would be of interest to extend these experiments into the localized regime, both 
in transmission and in reflection. This would require a substantial reduction in absorption, 
to ensure that L < £ < / a . 

Theoretically, much less is known about the delay-time distribution in transmission than 
in reflection. While we have a complete theory in reflection, as described in the previous sub- 
section, in transmission not even the N — 1 case has been solved completely. Regardless of 
the value of N, one would expect P(4>') for L ^> £ to have the same l/4>' 2 tail in transmission 
as it has in reflection, since in both cases the same resonances allow the wave to penetrate 
deeply into the localized region. For N = 1 this is borne out by numerical simulations by 



Bolton-Heaton et al. [pq] . These authors also used a picture of one- dimensional resonant 
transport through localized states to study the decay of the weighted delay time I(f>' (with 
/ the transmitted intensity). They found an algebraic decay for P(I<ft), just as for P((f>'), 
but with a different exponent —4/3 instead of —2. It is not known how this carries over to 
JV>1. 

Because of the finite length L of the waveguide, these algebraic tails are only an inter- 
mediate asymptotics. For N — 1 and exponentially large times |0'| > r s e L// ' the delay-time 



distribution has the more rapid decay |28|j2 



P(0')ocexp[-(//L)ln 2 (0'/r s )]. (17) 

Such a log-normal tail is likely to exist in the multi-mode case as well, but this has so far 
only been demonstrated in the diffusive regime I <C L <C £ |30 H32J1 . The l/<p' 2 intermediate 
asymptotics does not appear in that regime. 



III. HIGH-FREQUENCY DYNAMICS 
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FIG. 6. Distribution of the rescaled single-channel delay time 4>' = 4>' /4>', measured in trans- 
mission at a frequency v = oj/2tt = 18.1 GHz on the system described in Fig. 3. The smooth curve 
through the data is the analytical prediction (15) of diffusion theory (with Q = 0.31). From Ref. 

0- 



A. Reflection 

The high-frequency limit of the correlator @ of the reflection matrices is rather trivial. 
The two matrices r\oS) and r(u + Su) become uncorrelated for 5u> —>■ oo, so that C becomes 
the product of two independent random matrices. The distribution of each of these matrices 
may be regarded as uniform in the unitary group, and then C is also uniformly distributed. 
This is the circular ensemble of random- matrix theory [BO] , so called because the eigenvalues 
exp(i0 n ) are spread out uniformly along the unit circle. Their joint distribution is 



P({<f>n}) « II It* 

n<m 



— e 



4m 1/3 



(18) 



This distribution contains no dynamical information. 



B. Transmission 



The transmission problem is more interesting at high frequencies. Let us consider the 
ensemble-averaged correlator of the transmission matrix elements 



(a u (6u)) = (t* mn {uj)t mn {oj + 5u)). 



(19) 



10 



Following Ref. ||33| , we proceed as we did in Sec. [II A| , by mapping the dynamic problem 



without absorption onto a static problem with absorption. 

We make use of the analyticity of the transmission amplitude t mn (uj + iy), at complex 
frequency uj + iy with y > 0, and of the symmetry relation t mn {ui + iy) = t* mn (—uj + iy). 
The product of transmission amplitudes t mn {uj + z)t mn (— uj + z) is an analytic function of 
z in the upper half of the complex plane. If we take z real, equal to \5uj, we obtain the 
product t mn (oo+ |<5o;)tj! nn (u; — \5ui) in Eq. ( |T9D (the difference with t mn (uj + <5a>)£^ n (u;) being 
statistically irrelevant for 5uj <C uj). If we take z imaginary, equal to i/2r a , we obtain the 
transmission probability T = \t mn (uj + i/2r a )\ 2 at frequency uj and absorption rate l/r a . We 
conclude that the ensemble average of a can be obtained from the ensemble average of T by 
analytic continuation to imaginary absorption rate: 

M&*;)) = (T) for l/r a -► -iSu. (20) 

Higher moments of a are related to higher moments of T by (a p ) = (T p ) for l/r a — > —iSoj. 
Unfortunately, this is not sufficient to determine the entire probability distribution P(a), 
because moments of the form (a p a* q ) can not be obtained by analytic continuation. This is 
a complication of the transmission problem. The reflection problem is simpler, because the 
(approximate) unitarity of the reflection matrix r provides additional information on the 
distribution of the correlator of the reflection amplitudes. This explains why in Sec. [11 A 



we 



could use the mapping between the dynamic and absorbing problems to calculate the entire 
distribution function of the eigenvalues of r*(oj)r(oj + 5uj) in the limit 5u -^ 0. 

We will apply the mapping first to the single-mode case (N = 1) and then to the case 
iV ^> 1 of a mult i- mode waveguide. 

1. One mode 

The absorbing problem for N = 1 was solved by Freilikher, Pustilnik, and Yurkevich 
p3j |. Applying the mapping ( |2"0"| ) to their result we findQ 

(a w (&*;)) = exp(iSuL/c — L/l), (21) 

in the regime c/l <C 5uj <^C (cj 2 c//) 1 / 3 . (The high-frequency cutoff is due to the breakdown 
of the random-phase approximation fl35| .) The absolute value \(a)\ = exp(— L/l) is 5uj- 
independent in this regime. For L ^C I one has ballistic motion, hence (a) = exp(i6wL/c) is 
simply a phase factor, with the ballistic time of flight L/c. Comparing with Eq. (|2T|) we see 
that localization does not change the frequency dependence of the correlator for large Sou, 
which remains given by the ballistic time scale, but only introduces a frequency-independent 
weight factor. 



3 The coefficient in front of the factor L/l in the exponent in Eq. (21) would be — , according 
to the results of Ref. [[34|]. This would disagree with numerical simulations, which clearly indicate 
| (a) | = exp(— L/l) (K. J. H. van Bemmel, unpublished). The error can be traced back to Eq. (39) 
in Ref. @. 

11 



The implication of this result in the time domain is that (a w (t)) has a peak with weight 
exp(— L/l) at the ballistic time t = L/c. Such a ballistic peak is expected for the propagation 
of classical particles through a random medium, but it is surprising to find that it applies 
to wave dynamics as well. 

2. Many modes 

Something similar happens for iV 3> 1. The transmission probability in an absorbing 
multi-mode waveguide was calculated by Brouwer 



^ = N^kmrM -^i)- (22) 

for absorption lengths £ a = y/Dr & in the range I <C £ a <C £. The length L of the waveguide 
should be 3> I, but the relative magnitude of L and £ is arbitrary. Substitution of l/r a by 
— i8uj gives the correlator 

MSu)) = nlL17Xo. exp irm) ' (23) 



where td = L 2 / D is the diffusion time. The range of validity of Eq. (|23|) is L/£ <C \/r D 5u <C 
L/l, or equivalently D/^ 2 <C 5uj <C c/l. In the diffusive regime, for L -C £, the correlator 
(E3|) reduces to the known result [571 from perturbation theory. 



For max (D/L 2 , D/£ 2 ) <C <5co> <C c// the decay of the absolute value of the correlator is a 
stretched exponential, 

\M8w))\ = —yJr D 6wexp \-^~t d 8u - —J . (24) 

In the localized regime, when £ becomes smaller than L, the onset of this tail is pushed to 
higher frequencies, but it retains its functional form. The weight of the tail is reduced by 
a factor exp(—L/2Nl) in the presence of time-reversal symmetry. (There is no reduction 
factor if time-reversal symmetry is broken |33| . ) 

In Fig. 7 we compare the results of numerical simulations in a two-dimensional wave- 
guide geometry with the analytical high-frequency prediction. We see that the correlators 
for different values of L/£ converge for large 5uj to a curve that lies somewhat above the 
theoretical prediction. The offset is probably due to the fact that iV is not ^> 1 in the 
simulation. Regardless of this offset, the simulation confirms both analytical predictions: 



The stretched exponential decay ex exp(— JtdSuj/2) and the exponential suppression factor 
exp(— L/2£). We emphasize that the time constant To = L 2 /D of the high-frequency decay 
is the diffusion time for the entire length L of the waveguide — even though the localization 
length £ is up to a factor of 12 smaller than L. 

We can summarize these findings [Q for the single-mode and multi-mode waveguides 
by the statement that the correlator of the transmission amplitudes factorises in the high- 
frequency regime: (a^^Suj)) — > fi(Su>) /2(C) ■ The frequency dependence of f\ depends on 
the diffusive time through the waveguide, even if it is longer than the localization length. 
Localization has no effect on /1, but only on / 2 . 
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N = 5 



o L/l=18.1 
• L/l = 36.2 
□ L/l=?2.4 
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FIG. 7. Frequency dependence of the logarithm of the absolute value of the correlator (a^^duj)). 
The data points follow from a numerical simulation for N = 5, the solid curve is the analytical 
high-frequency result ( |24| ) for iV 3> 1. The decay of the correlator is given by the diffusive time 
constant Tp = L 2 /D even if the length L of the waveguide is greater than the localization length 
£ = 6 1. The offset of about 0.6 between the numerical and analytical results is probably a finite- N 
effect. From Ref. fp3l. 




IV. PROPAGATION OF A PULSE 

If the incident wave is a short pulse, then the separation into low- and high-frequency 
dynamics is less natural. Ideally one would like to know the entire time dependence of the 
correlator a^it) introduced in Eq. (p. A complete solution exists |38| for the ensemble- 
averaged correlator in the case of reflection, 



(au(t)) 



30 d5u 

c 



-iSujt 



«mM r mn(^ + <M) 



d5u> 



iSu)t 



N(N +1) J-oo 2n 
The second equality follows from the representation 

r(uj ± Su/2) = U T e ±i,s>/2 U. 



(TvC u {Su)). 



(25) 



(26) 



with $ = diag(0!, 2 , • • • 4>n) a diagonal matrix and U uniformly distributed in the unitary 
group. The factor 1 + 5 mn is due to coherent backscattering. It is convenient to work with 
the normalized power spectrum, 



N+l 1 r°° d5u 

1 + 6 mn u N J-oo 2ir 



iStot 



(TrC u (6u)), 



(27) 



normalized such that / °° dtV w (t) 
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Since e 1 ^" is an eigenvalue of the unitary matrix C, one can write 



2 77 



(Tr C u (6u))= / dtp^e*, (28) 



H 



where p(0) = (Z)n=i ^(0 — <f>n)) is the phase-shift density. This density can be obtained 
from the corresponding density p(R) of reflection eigenvalues R n (eigenvalues of rr^) in 
an absorbing medium, by analytic continuation to imaginary absorption rate: i/r a — > Su, 
R n — > exp(i(f> n ). The densities are related by 



/V 1 °° r l 

p(0) = — + -Re V e- in<t> / R n p{R) dR, (29) 

27T 7T „_-, JO 



as one can verify by equating moments. This is a quick and easy way to solve the problem, 
since the probability distribution of the reflection eigenvalues is known fL8|,19]| : it is given 



by the Laguerre ensemble ( jTT|) . The density p(R) can be obtained from that as a series 
of Laguerre polynomials, using methods from random- matrix theory [pOfl . Eq. (|29"D then 
directly gives the density p((f>). 

One might wonder whether one could generalize Eq. (EJ3) to reconstruct the entire distri- 
bution function P({0 n }) from the Laguerre ensemble of the i? n 's. The answer is no, unless 



5oj is infinitesimally small (as in Sec. [11 A| ). The reason that the method of analytic contin- 



uation can not be used to obtain correlations between the n 's is that averages of negative 
powers of exp(i<p n ) are not analytic in the reflection eigenvalues. For example, for the two- 
point correlation function one would need to know the average (exp(i0 n — i<f> m )) — > (R n R^) 
that diverges in the absorbing problem. It is possible to compute P({(p n }) for any 5u - 
but that requires a different approach, for which we refer to Ref. [|38| . 

The calculation of the power spectrum from Eqs. fl27|) — ( p9|) is easiest in the absence of 
time-reversal symmetry, because p(R) then has a particularly simple form. One obtains the 
power spectrum |38| 

V«{t) = -i^F(*/2JV 7 ), (30) 

1 "=*/t-l\ n „ ft 2 + l 



n=0 



m-— 1 E(— 1 )pJ— 1 , (3D 



where P n is a Legendre polynomial. (Recall that 7 = ar s , cf. Sec. [11 A| .) In the single-mode 
case Eq. (p0|) simplifies to 

K(t) = 2 7 (t + 2 7 )- 2 . (32) 

It decays as t~ 2 . For N — ► 00 Eq. ( pOf ) simplifies to 

^(t)=t- 1 exp(-t/ 7 )/ 1 (t/ 7 ), (33) 

where 1% is a modified Bessel function. The power spectrum now decays as t _3//2 . For any 



finite iV we find a crossover from V = J^/2-Kt 3 / 2 for r s < t C N 2 r s to V = 2N^t 2 for 
t ^> N 2 r s . This is illustrated in Fig. 8. 
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FIG. 8. Time dependence of the power spectrum of a reflected pulse in the absence of 
time-reversal symmetry, calculated from Eq. (BO) for N = 7 (open circles) and N = 21 (filled 



circles). The intermediate-time asymptote oc t -3 ' 2 and the large-time asymptote ex t~ 2 are shown 
as straight lines in this double-logarithmic plot. The prefactor is iV-independent for intermediate 
times but oc N for large times (notice the relative offset of the large-time asymptotes) . Courtesy 
of M. Titov. 



In the presence of time-reversal symmetry the exact expression for V^t) is more cumber- 
some but the asymptotics carries over with minor modifications. In particular, the large- N 
limit fl33|) with its t~ 3 ^ 2 decay remains the same, while the t~ 2 decay changes only in the 
prefactor: V = (N + l)^t~ 2 for t > N 2 t s . 

The quadratic tail of the time-dependent power spectrum of a pulse reflected from an 
infinitely long waveguide is the same as the quadratic tail of the delay-time distribution 
that we have encountered in Sec. O. It is natural to assume that the power spectrum for 
transmission through a localized waveguide of finite length has the same quadratic decay, 
with a cross-over to a log-normal tail for exponentially large times, cf. Sec. |LlCj 



V. CONCLUSION 



We now have a rather complete picture of the dynamics of a wave reflected by a disordered 
waveguide. The dynamical information is contained in the phase factors e l ^ n that are the 
eigenvalues of the product of reflection matrices r\u>)r(u> + 8u>). Three regimes can be 
distinguished, depending on the magnitude of the length scale ls u = J D/5u associated with 
the frequency difference 5u: 

• Ballistic regime, l$ u < I. This is the high-frequency regime. The statistics of the n 's 
is given by the circular ensemble, Eq. (|18|). 

• Localized regime, l^ > £. This is the low-frequency regime. The n 's are now dis- 
tributed according to the Laguerre ensemble (1121) . 

• Diffusive regime, I < Ig^ < £. The distribution of the n 's does not belong to any of 
the known ensembles of random-matrix theory |38 . 



15 



The emphasis in this review has been on the localized regime. The dynamics is then 
dominated by resonances that allow the wave to penetrate deep into the waveguide. Such 
resonances correspond to large delay times r n = lim^^o <p n /5uj. The distribution of the 
largest delay time r max follows from the distribution of the smallest eigenvalue in the Laguerre 
ensemble |[39|| . For /3 — 1 it is given by 



P(r max ) = jN(N + l)r m 4exp(-7iV(iV + l)/r max ). (34) 

It has a long-time tail oc l/V^ax? so that the mean delay time diverges (in the limit of an 
infinitely long waveguide). A subtle and unexpected consequence of the resonances is the 
appearance of a dynamic coherent backscattering effect in the distribution of the single-mode 
delay times. Unlike the conventional coherent backscattering effect in the static intensity, 
the dynamic effect requires localization for its existence. The recent progress in time-resolved 
measurements of light scattering from random media, reported at this meeting M, should 
enable observation of this effect. 

Extension of the theory to two- and three-dimensional localization remains a challenging 
problem for future research. We believe that the dynamic coherent backscattering effect 
will persist in higher dimensions, provided the localization length remains large compared 
to the mean free path. Several methods have been proposed to distinguish absorption 
from localization in the static intensity p0| , |4l| . The effect reviewed here could provide this 



information from a different, dynamic, perspective. 
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